Aspect ratio dependence of charge transport in turbulent electroconvection 
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We present measurements of the normalized charge transport or Nusseh number Nu as a function 
of the aspect ratio F for turbulent convection in an electrically driven film. In analogy with turbulent 
Rayleigh-Benard convection, we develop the relevant theoretical framework in which we discuss 
the local-power-law-scaling of Nu with a dimensionless electrical forcing parameter TZ. For these 
experiments where lO" < 7^ < 2 x 10^ we find that Nu ~ F(r)7l^ with 7 either = 0.26 (±0.02) 
or 0.20 (±0.03), in excellent agreement with the theoretical predictions of 7 = 1/4 and 1/5. Our 
measurements of the aspect ratio-dependence of Nu for 0.3 < F < 17 compares favorably with the 
function F(F) from the scaling theory. 



Heat transport has for many years been the corner- 
stone in the study of turbulent Rayleigh-Benard convec- 
tion (RBC) IJ. Over the past three years there has been 
remarkable experimental [3) S S IE 0. 13 and theoret- 



ical 0, Ei progress in characterizing the properties 
and mechanisms of heat transfer in fluids heated from 
below. So far, experiments and theory are relevant only 
to approximately unit aspect ratio systems i.e. geome- 
tries that have comparable lateral and vertical dimen- 
sions. Turbulent convection flows in nature, however, are 
laterally extended and the transport of heat across a layer 
of fluid has been inadequately studied as a function of the 
aspect ratio F = lateral dimension/ vertical dimension. 
In particular, experiments at F > 2 have seldom been 
performed primarily because it is extremely difficult to 
achieve the strong forcing that is readily reached in F « 1 
containers. 

In this Letter, we present a study of an analogous elec- 
trically driven convecting system which allows a broad 
range of F to be explored at moderate levels of forcing. 
We show that measurements of the normalized charge 
transport Nu varies with the aspect ratio F consistent 
with a function F given by the scaling theory. The strong 
dependence on F especially for F « 1 implies that exper- 
iment and theories relevant to this regime are restrictive. 
The data and the function F become increasingly inde- 
pendent of F for F ^ 1. This observation highlights the 
importance of developing experiments and theory for lat- 
erally extended systems where it appears that universal 
behavior, often lacking in F « 1 systems, may be re- 
stored [131 ■ Our data also reveal that Nu ~ TZ'^ with 
7 = 0.26 and 0.20 in agreement with a local-power-law- 
scaling theory developed in a manner identical to the 
Grossmann-Lohse (GL) model for turbulent RBC 

In previous work, it has been demonstrated that an 
electrically conducting, freely suspended liquid crystal 
fllm between parallel wires could be driven to convect 
when a sufficiently lar ge potential drop was ap plied 
across its edges 
periments ll 
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FIG. 1: Schematics of the experiment: (a) top and (b) side. 

extended to the naturally periodic geometry of an annu- 
lar film shown schematically in Fig.^ Analogous to RBC 
where an inverted mass density distribution is unstable 
to buoyant forcing, the thin film has an inverted surface 
charge density distribution that is unstable to electrical 
forcing. In this system, an initial bifurcation to convec- 
tion rolls occurs when the applied voltage V exceeds a 
critical voltage Vc, corresponding to the critical temper- 
ature difference in RBC. Secondary bifurcations occur 
at higher forcing and result in changing the number of 
counter-rotating vortex pairs that are arranged around 
the annulus. At much higher forcing the fiuid becomes 
turbulent but retains the large scale structure of the rolls 
or convection cells. 

Electroconvecting smectic films have several advan- 
tages and disadvantages relative to traditional RBC ex- 
periments. In turbulent RBC experiments, the accurate 
measurement of heat transport requires detailed account- 
ing of the heat conducted through the sidewalls 3. The 
annular geometry of the smectic film (see Fig. ^) is 
free of lateral boundaries and consequently without side- 
wall losses, which facilitates precise measurement of the 
charge transport between the inner and outer electrodes. 
The film is an annular disk of width d — — ri ~ 5 
mm and thickness s ~ 0.1 /xm. Here (tq) are the 
radii of the inner (outer) electrodes that support the film 
(see Fig. ^) . These radii can be varied so as to achieve 
different values of F. The tiny size of the film, which con- 
tains many orders of magnitude less working fiuid than 
required for an RBC experiment, results in data acqui- 



sition timescales of minutes, rather than days. On the 
other hand, the smectic film is dehcate and while many 
films survive vigorous forcing without thickness change, 
some do suffer sudden thickness variations. Since Vc and 
other properties of the flows depend on thickness, these 
data must be discarded. The dc electrical forcing also 
results in conductivity drifts in the films, which result in 
systematic uncertainties. The experiments we describe 
here were carried out at atmospheric pressure, which in- 
creases the film stability. But air drag, which is known to 
produce a quantitative but not qualitative change to the 
dissipation, is not properly accounted for in the theory. 
Nevertheless, notwithstanding these caveats, we are able 
to efficiently explore the turbulent scaling regime over a 
broad range of F, complementing RBC experiments. 

Our experiment consists of a temperature controlled 
film of octylcyanobiphenyl (8CB), a smectic- A liquid 
crystal, suspended between two concentric gold-plated 
electrodes. In the smectic-A phase, the film flows as 
an isotropic, incompressible and Newtonian fluid in the 
plane of the film. The flow is strictly two-dimensional. 
The film is driven to electroconvect by a dc voltage ap- 
plied to its inner edge while holding the outer edge at 
ground potential. The annular assembly is enclosed in 
a Faraday cage. An experiment consists of imposing a 
voltage V to the inner electrode and measuring the cur- 
rent / transported through the film with a sensitive elec- 
trometer. The applied voltage was varied between and 
1000 volts in a sequence of small incremental and decre- 
mental steps resulting in a current-voltage characteris- 
tic. Further details of the apparatus and procedure can 
be found in Refs. [H El El EH. Figure M shows 
a representative current-voltage (IV) characteristic. For 

V < Vc the fluid is quiescent with the current sustained 
by ohmic conduction. When V > Vc the fluid is orga- 
nized in counter-rotating vortices and additional charge 
is transferred by convection; this is seen by the increase 
in slope of the IV-characteristic. At higher voltages, the 
fluid becomes turbulent and the transition is marked by 
a sudden increase in the rms fluctuations of the current, 
as shown in the inset of Fig. 13 

A film is a 2D annular sheet. The film has thick- 
ness/width s/d ~ lO^**. Geometrically, a film is de- 
scribed by its radius ratio a — ri/ro and the aspect ra- 
tio Tr = 7r(ri + ro)/{ro — ri), which is the ratio of the 
mid-radius circumference to the film width. So defined, 
TT < Fr < oo. To make correspondence with the con- 
ventional aspect ratio V for RBC, we define F = F^ — tt 
for the annulus. Two other dimensionless parameters 
describe the experimental system: TZ = e^V^ /ar]s^ and 

V = eorj/pasd. Here TZ is the control parameter and is 
a measure of the external electrical forcing and V, the 
Prandtl-like parameter, is the ratio of the time scales of 
electrical and viscous dissipation processes in the film. 
In the above the fluid density, molecular viscosity, and 
conductivity are denoted by p, rj, and cr, and eg is the 
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FIG. 2: A representative current-voltage curve for an annu- 
lar film showing the onset of electroconvection. Data obtained 
for increasing (decreasing) voltages are shown in A(v). The 
inset shows the rms fluctuations of the current vs. the appfied 
voltage. The transition to turbulent flow is marked by a sud- 
den large increase in the fluctuations. Here F = 6.60 ± 0.02 
and P = 36 ± 1. 

permittivity of free space. 

We normalize the measured total electric current / by 
the portion due to conduction Icond = cV where c is 
the ohmic conductance of the film obtained from the IV- 
characteristic when the film is not convecting, i.e. when 
V < Vc. We define the normalized charge transport or 
Nusselt number Nu = I /Icond- Figure 13 shows repre- 
sentative sets of Nu vs TZ data at F = 0.33 ± O.OI and 
6.60 ±0.02. We find that Nu - W for 10^ <TZ<2x 10^ 
with either 7 = 0.26 ± 0.02 or 7 = 0.20 ± 0.03, depend- 
ing on the size of V. The error bar is obtained from the 
variation in the best-fit value of 7 over 26 data sets for 
all four values of F. 

We developed a theoretical model describing turbulent 
electroconvection by borrowing from turbulent RBC. The 
equations of motion that describe electroconvection con- 
sist of the incompressible Navier-Stokes equation supple- 
mented with an electrical body force, the charge conser- 
vation equation that accounts for ohmic conduction and 
advection of charge and Maxwell's equation that relates 
the surface charge density to the electric potential; 



i9tu + u • Vu = 
dtq + u-Vq = 

q = 



-V- + t/V^u - -VV', 

p p 

-26092-0 . 



(1) 

(2) 
(3) 



These equations are constrained by the no-slip and ap- 
plied electric potential boundary conditions. In the above 
equations V, u, p and v are the 2D gradient, velocity, 
pressure and kinematic viscosity respectively. The elec- 
tric potential -0 is three dimensional and extends outside 
the film, where (V^ -I- dDip = 0. The Maxwell equation 
relates the perpendicular gradient of ip to the surface 



3 



charge q. The factor of 2 arises from the film's two free 
surfaces. The relation between ip and q is thus nonlocal 
and somewhat complex. It can be considerably simpli- 
fied by making a local approximation, setting q = 2/3'0, 
where /3 is a certain constant. This local approximation 
has been shown to be adequate in describing the onset of 
electroconvection. See Ref. |0 for a detailed discussion 
of the theoretical model. 

The Eqns.^lSare similar to the Boussinesq equations 
for turbulent RBC. In particular, in the local approxima- 
tion, there is an almost precise correspondence between 
the temperature and the electric potential. The period- 
icity of the annular geometry and thereby of the velocity 
and electric potential allow for exact relations for the 
globally-averaged kinetic and electric dissipations. De- 
noting averages over the system volume by (• • •) we find 
the following for the kinetic dissipation eu = (i^(Vu)^} 
and for the electric dissipation = (ct(VV')^): 



ln(l/a)(r2 - rf){ro ~ 
ln(l/a)(r2 — r^) 



(4) 
(5) 



These relations are similar to those for the kinetic and 
thermal dissipation in RBC presented in Refs. 9. 10. 22,]. 
Following the GL theory, we decompose the kinetic and 
electric dissipations into bulk and boundary layer parts: 
eu.i/; = Cu V"^^ + V"- Considering the various combina- 
tions of the dominant contributors to the total dissipa- 
tion (e.g. Tilj-ip and BULK-u etc.) we can determine 
the corresponding local power-law scalings for each of 
the regimes. We assume that the turbulent flow is com- 
prised of convection cells that tile the annulus as shown 
in Fig. ^ and that each cell is roughly square i. e. it has 
the same transverse and lateral dimension. The large 
scale circulation is then the vortex that defines a cell. 
Assuming laminar boundary layer scaling we can then 
show that 



Nu - F(^)7^^7'^ Re-7^''*P**, and (6) 

(7) 



F(r) ^ £±^ln 



As might be expected, we find the same set of expo- 
nents 7,(5, 7*, (5, that appear in the GL theory. Here, Re 
is the Reynolds number of the large scale circulation. 

In the relatively small TZ regime, the total dissipation is 
dominated by the BL-contributions and the theory pre- 
dicts 7 = 1/4 and 6 = 1/8. In a neighboring regime 
where the dissipation is primarily e^^'"^ and e^^, the 
theory predicts j = S = 1/5. Our measured 7 exponents, 
as shown in Figure|2|are in reasonable agreement with the 
theoretical predictions from the GL model, if we suppose 
that we traverse the appropriate regimes as V varies. A 
detailed validation of the theory would depend crucially 




FIG. 3: Plots of Nu vs 7^ for F = 0.33 ± 0.01(A) and 6.60 ± 
0.02(0). For these V = 8.8± 0.3 and 36 ± 1 respectively. The 
solid reference lines have slopes of 1/5 and 1/4 respectively. 

on either directly demonstrating the correctness of the 
assumptions about the dominant contributions to the dis- 
sipation or on systematically showing that the right local 
power-law scalings are indeed obtained as one varies V 
and TZ. Our experiments currently span only relatively 
low TZ < 2x 10^ and sparsely cover the rather wide range 
5 < < 250. Many more experiments will be needed to 
adequately test the GL theory for turbulent electrocon- 
vection. 

In the theoretical treatment described above, we have 
taken into account the dependence on the aspect ratio 
of the system, a consideration that was missin g in the 
GL theory which treated only F « 1 systems [sl, lnJ, |llj . 
We find that the charge transport is modified by a F- 
dependent factor F(F). Since the large scale circulation 
is local to each convection vortex. Re is independent of 
F. Unlike Ref. 22], we find the aspect-ratio dependence 
is not a power-law scaling but rather a function of the fi- 
nite annular geometry. We can make a direct comparison 
with previous turbulent RBC experiments by making a 
correspondence between the different cell geometries, us- 
ing F = F,. — TT as described above. We plot fcF(F) vs. 
F, with the constant k = n/{{n + 1) ln(27r + 1)) = 0.382. 
Following Ref. 3]j we choose this normalization so that 
A:F(F = 1) = 1. Then kF 0.764 in the F ^ 00 limit. 
The function fcF decreases monotonically with F with 
the greatest variation for F < 2, and is within 2% of the 
limiting value for F > 7. 

Our experimental data span the range 0.3 < L < 17. 
From power-law fits to Nu vs. TZ data we have extracted 
the exponents 7. To determine F(F) we then divide 
Nu by 7^''7"^. Because our data extend over a rather 
wide range of we expect to cross regimes with differ- 
ing 7 and (50, 0]. To extract the aspect ratio depen- 
dence alone, we used the fitted value of 7, and where 
7 w 0.25 (0.20), we used the GL-theory prediction of 
(5 = 1/8 (1/5). Using one free parameter for all the data. 
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FIG. 4: The function kF vs F. fcF(F) is normalized to unity at 
F = 1. The solid circles (•) are from the present work. Data 
from RBC experiments in AcetonejWith (O) and without (□) 
sidewall correction are from Refs. Q and Q]. Also shown are 
data in Air (A), from Ref. [H, Helium (v) from Ref. ^ 
and (<) from Refs. 1,'^ and water {>) from Refs. The 
inset shows the same data on a logarithmic scale. 



we again scaled these results so that F(r = 1) = 1. Our 
data for 4 different T are in reasonable quantitative agree- 
ment with the theoretical function F, as shown in Fig.^] 
Each data point is an average over 4—10 runs. The error 
estimates are representative of the scatter in F over these 
data. In particular, the largest contribution to the error 
at r = 16 is systematic and arises from the ambiguity in 
the correct exponent for V. 

Data from several turbulent RBC experiments0,0,0j 
1^ for 4 values of F are also broadly in agreement 
with the function F, in spite of the difference in geometry 
and the higher range of Rayleigh numbers. Earlier RBC 
experiments |25j used gases as the working fluid and 
deviate significantly from F. 

Several important questions are raised by our results 
on the aspect ratio dependence of charge transport in tur- 
bulent electroconvection and, by extension, that of heat 
transport in turbulent RBC. The strong dependence of 
F on F for < r < 7, suggests that the dynamics of 
turbulent convection is dependent on the system's lat- 
eral extent. That the charge or heat transport, a global 
property, is sensitive to the aspect ratio further suggests 
that local properties such as temperature and charge fluc- 
tuations may be even more strongly dependent on F. Ex- 
periments in turbulent RBC clearly show that while heat 
transport is insensitive to cell geometry at fixed F « 1, 
the fluctuations are strongly affected by the shape of the 
lateral boundary ■ The relative independence of F on 
F for large aspect ratio, F > 7 suggests that the nor- 
malized charge and heat transport approach a universal 
value in laterally extended systems. We conjecture that 
fluctuations in the interior may also become universal in 
large aspect ratio systems. 



Our results emphasize the importance of extending 
convective turbulence experiments to larger aspect ra- 
tios. For turbulent smectic electroconvection, we would 
like to make a systematic study of the scaling of Nu in 
the parameter space of TZ and V. V can be varied in 
principle by changing the thickness and/or width of the 
film. We would also like to study long time series of Nu in 
order to infer the characteristic size and duration of the 
fluctuations. It is a challenge to extend the forcing pa- 
rameter TZ to values much beyond 10^, since this requires 
applying large voltages across the film that may result in 
dielectric breakdown. A first step would be to better 
control the electrochemistry of the liquid crystal film to 
reduce drifts and push down the critical voltage Vc- It 
would also be interesting to extend oiir previous studies 
of electroconvection under shear 0, 0, 0, |^ to 
the turbulent regime. Finally, we would like to develop 
local probes to study the fluctuations of the velocity and 
electric potential over the film. 
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